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Abstract 

We describe an experiment involving a mass oscillating in a viscous fluid and analyze viscous 
damping of harmonic motion. The mechanical oscillator is tracked using a simple webcam and an 
image processing algorithm records the position of the geometrical center as a function of time. 
Interesting information can be extracted from the displacement-time graphs, in particular for the 
underdamped case. For example, we use these oscillations to determine the viscosity of the fluid. 
Our mean value of 1.08 ±0.07 mPas for distilled water is in good agreement with the accepted 
value at 20° C. This experiment has been successfully employed in the freshman lab setting. 
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I. INTRODUCTION 



The topics of oscillatory and simple harmonic motion (SHM) are of fundamental im- 
portance to physics and engineering. The differential equations that lie at the heart of 
oscillatory motion are frequently covered in introductory undergraduate courses. The appli- 
cations are also widespread, touching upon phenomena as diverse as vibrations of atoms in a 
crystal, current in LCR circuits, the transmission of electromagnetic radiation in dielectrics 
and chaos. However, we felt that even though a large number of our students knew how 
to solve these second order differential equations, very few had actually seen the solutions 
naturally emerge in real physical situations, especially in the formal laboratory environ- 
ment. The present experiment demonstrates these solutions and provides physical insight 
into light and heavy damping. The data, that is acquired from an off-the-shelf webcam, is 
subsequently processed to calculate a useful parameter, the viscosity of commonly available 
fluids such as distilled water, ethanol and methanol. The values are in reasonable agreement 
with accepted values. 

In the process, students are also exposed to techniques in image processing. The use of 
these techniques is becoming more widespread especially in the area of video microscopy 
where a CCD camera is used to track the motion of particles, microspheres or fluorescent 
proteins in fluidic environments. Some of these ideas have also been discussed^ in this 
journal and used to follow Brownian motion and its dependence, for example, on fluid vis- 
cosity, particle size and temperature as well as the estimation of Boltzmann's and Avogadro's 
constants^. The present experiment may serve as a first introduction to some commonly used 
algorithms and tools in image processing such as frame grabbing, color control, motion track- 
ing and using videos for the quantitative verification of mathematical and semi-empirical 
models. 

Not surprisingly, there is already a diverse repertoire of experiments performed in in- 
structional physics laboratories worldwide that analyze different facets of harmonic motion. 
These include the use of oscillating water columns^, swinging pendulums^ and masses at- 
tached to springs. These experiments analyze SHM as well as its nonlinear generalizations 
using potentiometers^, photocells^, photosensors^ and force sensors^. In the present exper- 
iment, we use a webcam to track the damping of a mass oscillating in various fluids such 
as honey, water, ethanol and methanol and adjust the fluids to exhibit heavy and light 
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damping. 

Several illuminating articles revolving around the concept of viscous damping have ap- 
peared from time to time in pedagogical journals. In our case, we process the images 
acquired from the webcam, in the case of light damping, and go on to calculate the absolute 
viscosity of the fluid, a parameter that is conventionally determined using viscometers of 
various kinds 9 . Our proposed technique is a new method of acquiring data from a viscously 
damped system and differs from traditional methods in several useful respects. For example, 
Alexander and Indelicate^ used a force sensor to monitor the damping of a spherical mass 
in water and processed their data to determine the viscosity. Their calculation differed from 
the accepted value by an order of magnitude, indicating the presence of large, systematic 
errors. Another article^ described the use of a photosensor to track the viscous damping of 
a pendulum in air. This method works very well for gases but the downside is that ordinary 
photosensors cannot be immersed in liquids. Our method employing the webcam acquires 
data remotely obviating the need for any physical contact with the oscillator. The posi- 
tions are directly measured and the time stamped displacements can also be numerically 
processed, if desired, to estimate the velocities^. We show that in conjunction with the 
appropriate mathematical model, ours is a quantitatively accurate method for determining 
the viscosity of liquids. 

The present experiment extends the list of webcam-based experiments oriented for un- 
dergraduate labs demonstrating the diffraction of light, the diffusion of ink in water 11 , as 
well as quantitative measurements on shadows, sprouting water jets, hanging chains and 
caustic reflections^. 

II. THEORY 

The physical model is of a linear harmonic oscillator, a spherical mass attached to a spring 
that remains within its elastic limit and executes simple harmonic motion inside a viscous 
medium that we call the fluid. At small speeds, a good approximation to the damping force 
is given by^, 

F « bv, (1) 

where v = y is the linear speed of the oscillator and b is a proportionality constant charac- 
terizing the medium and the shape of the oscillator. The characteristic equation of motion 



3 



is. 

,2 



y + "/y + uy = £, (2) 

where 7 = b/m, u 2 = k/m, k is the spring constant, m is the mass of the particle and £ is 
the ratio of the buoyant force to the mass of the sphere. The constants 7 and u represent, 
respectively, the damping constants and the natural frequency of the oscillator when the 
damping is switched off. During the experiment, the oscillator is completely immersed in 
the fluid at all times. As a result, the buoyant force is constant and always directed upwards. 
Therefore, £ is a constant and causes a uniform offset in the displacement y; mathematically, 
this is a shifting of the solutions of Equation (J2J) on the vertical axis (in a displacement-time 
graph) . 

A general expression for the damping force experienced by a spherical mass of radius r, 
oscillating in a fluid of viscosity rj and density p is given by^^i, 

(t \ / 2 t \ dv 

1 + _j„_ 3 ^l + 5 -j M _, (3) 

where 5 is called the penetration depth, 

S - M. (4) 
V 

The penetration depth defines the depth at which the wave amplitude falls to 1/e (about 
37%) of the amplitude at the surface of the sphere. This depth decreases with the frequency 
of the wave, but increases with the kinematic viscosity rj/p of the fluid. From Equation Q, 
we recover Stoke's famous law, F = —Qnrjrv, if the mass is falling with uniform velocity and 
uj is zero. If the dimensions of the sphere are greater than the skin depth, r^>5, Equation (J3]) 
reduces to, 

_, 67r??r 2 2irr 3 dv , , 

' v-— p-771 (5) 
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but the expression applies if the mass is oscillating in a container of infinite dimensions; for 
the physically realizable case of a finite container, a correction term is added^, 

/ 67T7]r 2 2irr 3 dv\ f r 
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d being the radius of the cylindrical container. Now from ([!]) and ([6]), it can be shown that 



the oscillating mass obeys the equation of motion fl2]) with the mass m replaced by m f and 
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From Equation (jSJ) we deduce an interesting fact: if the container is infinitely large, exactly 
half of the fluid displaced by the sphere accompanies it in its oscillatory journey by forming 
a boundary film around the sphere! We estimate the damping coefficient from the mass that 
negotiates underdamped SHP and use the relationship ((7j) to estimate the viscosity. Finally, 
as a simple demonstration, students employ solutions, such as honey in water, adjust their 
concentrations and investigate light and heavy damping. 

III. THE EXPERIMENT 

The experimental setup is simple, easily repeatable and is illustrated in Figured] A sphere 
of mass 210 g is attached to a spring of spring constant 8 Nm" 1 and set to oscillate vertically 
in the fluid while a webcam acquires a time-series of frames. In our case, the nominal frame 
acquisition rate is set at 30 per second with 300 frames acquired in total. A white sheet of 
cardboard is placed in the background to help reduce noise levels while the spherical mass is 
also painted black also that helps in further reducing the light reflections. Figure EJa) shows 
a string of acquired colored frames. To reduce the noise levels even further, the frames are 
cropped showing only the interesting region in which the mass oscillates. The image and 
data processing are all done in Matlab. However, this is only out of convenience and the 
students' prior training — any other commercially available software can be used with equal 
ease. The data analysis is described in the next section. 

After the frames have been acquired, the next step is to calculate the centroid of the 
oscillating sphere. This a common technique in image processing and Matlab provides an 
algorithm for the purpose 17 . Since the accurate determination of the centroid is susceptible 
to color fluctuations, the frames are converted to black and white at this stage with the 
mass appearing white against a black background. Figure EJ^b) shows a string of black and 
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FIG. 1: Schematic diagram of the experimental set-up. The apparatus is adjusted so that the 
mean position of the spherical mass is in the horizontal line of sight of the webcam. 

white, cropped frames, corresponding to the colored ones in Figure [2(a). The centroids of 
the white regions in Figure [2(b) are determined for all the 300 frames and then plotted as 
a function of time as we now discuss. 



FIG. 2: (a) Colored photographs of the frames showing the oscillations, (b) Typical black and 
white frames after processing. The images are cropped to show only the region of interest. 



IV. DATA ANALYSIS 

The solutions of the differential equation ([2]) reveal three different kinds of harmonic 
motion depending on the relative magnitudes of 7 and u. For example, the condition u> 2 < 
7 2 /4 represents overdamping (or heavy damping). We experimentally achieve this regime by 
oscillating the mass in a jar of honey whose density is empirically adjusted between 1.28 and 
1.35 g cm -3 . This case is exhibited in Figure [3j showing the displacement y as a function 
of time, the curves (a) through (c) represent increased concentrations and densities of the 
honey mixture, resulting in heavier damping. The viscous drag of the thicker fluid causes the 
mass to only slowly approach equilibrium, without performing any oscillations. As the plots 
show, the approach towards equilibrium becomes slower from (a) to (c), i.e., with increasing 
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damping coefficient. The general form of the heavily damped solution is given by^, 

y(t) = Ae-^ + Be~^\ (10) 

where 

^1,2 = ^7T ^7 2 -^ 2 . (11) 

The constants A and B have dimensions of length and depend on the initial state (position 
and velocity) of the mass, where t = is defined by the registration of the first frame. In 
our case, A and B are determined from fitting the displacement-time curve to the solution 

(HUD- 
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FIG. 3: Damping in a jar of honey depicting overdamped SHM and the varying rates of approach 
towards equilibrium. The parameters are (a) m = 1.475 s _1 , \ii = 0.002 s _1 , A = 69 p and 
B = 39 p; (b) m = 0.711 s -1 , \i % = 0.021 s~\ A = 54 p and B = 51 p; and (c) m = 0.484 s -1 , 
H2 = 0.025 s _1 , A = 59 p and B = 64 p. 

The case of light or underdamped SHM, 7 u; is achieved when the mass is made to 
oscillate in a less viscous medium such as distilled water or pure ethanol. The angular 
frequency of the mass cu w is then given 



7 



FIG. 4: (a) Damped simple harmonic oscillations versus time. The fit parameters are A = 
21 A pixels, 4> = 3.82 radians and the decay constant, 7=0.1398 s~~ . (b) Fluctuations in area 
of the white sphere versus time. 

and the corresponding displacement becomes, 

y(t) = Ae~ lt/2 cos(co w t + 0), (13) 

where A and <fi are, respectively, the amplitude and phase determined from the initial con- 
ditions. 

Figure H^a) shows the experimental data points along with the predicted theoretical 
model, encapsulated by the expression ffTB"]) . The data points are fit to the solution and 
the parameter 7, representing the decay of the oscillations, is determined. Finally, the 
viscosity is calculated from Equation ([7]). For distilled water, our value of rj equals (1.08 ± 
0.07) mPas and is in good agreement with the accepted value of 1.002 mPas at 20°C 
For ethanol and methanol, the determined viscosities are (1.20 ± 0.07) mPas and (0.39 ± 
0.07) mPas, concordant with the viscosities determined by other methods^ 1 ^, 1.1005 mPas 
and 0.53 mPas. 

The uncertainty in the results can principally arise from two sources: uncertainties in 
time and in position. With a digital camera, frames are dropped or delayed and some 
frames can even arrive earlier than expected. These variabilities give rise to timing jitter. 
We circumvent this problem by using the real time stamps instead of the nominal inter- 
frame delay of 1/30 s (inverse of the frame rate)^. In case a frame is dropped or delayed, 
the time data explicitly shows that. The time stamps for a typical run of the experiment 
are shown in Figure and reveal the inter-frame variation. Finally, the least count of the 
time measurements as resolved by our software is ±1 millisecond. 

The uncertainty in position can be caused by multiple sources. One source is the spurious 
reflections of light from the metallic sphere, the glass of the beaker, refraction and the 
lighting conditions inside the room. For example, while computing the centroids in the 
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FIG. 5: Plot of the time lapse between successive frames in a typical run of the experiment. It is 
clear that the recorded frames cluster in two groups. One is the group of frames that are acquired 
at or near the nominal frame rate of 1/30 ~ 0.033 s whereas the second is the group of frames that 
are delayed by twice the sampling time, ~ 0.066 s; the latter represent dropped frames. We use 
the actual time values in our experiment instead of an average time. 

acquired frames, it was noticed that the 'white' area fluctuates due to scattering. These 
aerial fluctuations are shown in Figure 5(b), with a scatter of ~ 5 A = ±45 pixels 2 . For an 
average radius r ~ 32 pixels and an average area A = irr 2 , the uncertainty in the radius is 
Sr = r5A/{2A) « 0.2 pixels. 

In addition to the vertical motion, the oscillator also swings horizontally, making the 
system resemble a pendulum with periodically oscillating length, a so-called parametric 
oscillator 21 . In Figure H](a), we depict the horizontal displacement superposed on top of the 
vertical oscillations. The motion for such a system deviates from simple harmonic and in fact 
becomes nonlinear, governed by nonlinear differential equations (e.g., Equations 1,4,6 ir*22) 
that are solved numerically. Then our predicted solution (1131) should only be approximate. 
However, the excellent curve fit indicates that the nonlinearities are indeed small and the 
damped oscillations represent, to a fairly good extent, en exponentially damped sinusoid. 

The horizontal amplitudes are small when compared with the vertical ones, approximately 
one part in ten (10%). We can use these small horizontal amplitudes, 5x, to estimate the 
uncertainty in the vertical coordinates, 5y, under certain simplifying assumptions. Suppose 
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9 is the angle subtended by the horizontal swing of the pendulum from the vertical. If Iq 
is the mean length of the spring-mass system and Iq ± Al max are the length extrema, the 
angles at the extrema are given by 9 = x/(Iq ± Al max ) and for small 9, we have x ~ x, 
9 w x/{Iq ± Al max ). When both the spring extension and the horizontal swings are at 
their maximums, the centroid-determination algorithm records a vertical position y = (l + 
Al max ) cos9 whereas the position in the absence of horizontal motion would be Iq + Al max , 
the deviation being 5y = (l Q + Al max ) (cos 9 - 1) = (/ + A/ ma:E )(2 sin 2 6»/2) « (/ + A/)(6» 2 /2). 
With the small angle approximation, 

^ ~ / +A7 ' (14) 
and similarly, for the maximum compression, 

5y ~T^Xl — ' (15) 

With ball-park figures of Iq ~ 200, A/ maa; rj 20 and 5x ~ 2 (all units in pixels), the 
uncertainty Sy is < 0.02 pixels. In practice, the horizontal and vertical displacements are 
not in phase and the actual uncertainties will lie in the range specified by (114j) and (1151) . 
Therefore, based on the light scattering and horizontal motion, a conservative guess on the 
order of magnitude of the positional uncertainty would be ±1 pixels. 



V. CONCLUSIONS 



In this experiment, we have used a new technique of employing a webcam to acquire 
time stamped frames of a mass oscillating inside a viscous fluid. Our students take keen 
interest in employing the webcam, generally a household object of entertainment value, for 
a scientific pursuit. The experiments are low cost, ideal for cost-conscious experimental labs 
and easily reproducible. Furthermore, the experimentally determined viscosities of water, 
ethanol and methanol are in excellent agreement with previously reported values. Notably, 
our method returns a value far more accurate than a similar experiment that employed a 
force sensor-. 

The experiment introduces students to damping and simple harmonic motion by acquiring 
data remotely. Students learn to work in pixel co-ordinates, perform simple image processing 
and computational tasks and test experimental observations with quantitative predictions, 
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in our case, solutions of second order differential equations. One can also use our webcam- 
based approach to study the phenomenon of parametric resonance 23 . Thus we can oscillate 
the combined mass-spring and pendulum system with different excitation frequencies and 
study the (in) stability properties of the nonlinear system as well as its transition into chaos. 
This is an interesting idea to pursue for short-term research or a lab project. 

The experiment has been performed by a freshman class of about 150 students and very 
well received. It is hoped that the new method of analyzing simple harmonic motion will 
be welcomed by physics teachers. The authors like to thank Waqas Mahmood, Muhammad 
Wasif and Umer Suleman for demonstrating the experiment to the class. 
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